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ABSTRACT 

Observations of clusters of galaxies show ubiquitous presence of X-ray cavities, pre- 
sumably blown by the AGN jets. We consider magnetic field structures of these cav- 
ities. Stability requires that they contain both toroidal and poloidal magnetic fields, 
while realistic configurations should have vanishing magnetic field on the boundary. 
For axisymmetric configurations embedded in unmagnetized plasma, the continuity 
of poloidal and toroidal magnetic field components on the surface of the bubble then 
requires solving the elliptical Grad-Shafranov equation with both Dirichlet and Neu- 
mann boundary conditions. This leads to a double eigenvalue problem, relating the 
pressure gradients and the toroidal magnetic field to the radius of the bubble. We 
have found fully analytical stable solutions. This result is confirmed by numerical 
simulation. We present synthetic X-ray images and synchrotron emission profiles and 
evaluate the rotation measure for radiation traversing the bubble. 

Key words: MHD - methods: analytical - methods: numerical - galaxies: clusters: 
intracluster medium - plasmas, X-rays: galaxies: clusters. 



1 INTRODUCTION 

X-ray and radio observations of intracluster medium (ICM) 
demonstrate the presence of cavities produced by the AGN 



jets flMcNamara fe Nulsen||2007| |Fabian et aL|2000| |Birzan 



|et al.||2004| |Diehl et al.||2008| |Dong et al.||2010| ). Cavities 



are filled with jet-supplied plasma and magnetic fields and 
are observed both as depressions in X-ray surface bright- 
ness and through intrinsic synchrotron emission. The na- 
ture of the plasma inside the cavity is not determined yet, 
a way to determine whether the plasma is thermal or not is 

|2005) . 



through Sunyaev-Zel'dovich effect JPfrommer et al. 



Most cavities are not expanding supersonically and, thus, 
are in pressure balance with their surroundings. Their rise 
in the cluster potential is caused by buoyancy. 

Cavities are long lived entities, surviving on time scales 
much longer than the time scale of buoyant rise. On the 
other hand, fluid (non-magnetic) simulations show that light 
bubbles are disrupted on one scale heig ht (|Robinson et al. 
2004] IChurazov et al.||200T] [Bruggen fe Kaiser ||2002| | Jones 



&; De Young 2005), but there have also been works which 
suggest that purely hydrodynamics effects stabilise the cav- 
ities ( |Reynolds et al.||2005| |Scannapieco & Bruggen 2 010| ). 
Magnetic field can stabilize the cavities via two mechanisms. 
First, magnetic field in the ICM drapes around expanding 
cavities ( Lyutikov|2006 ) forming a layer of nearly equiparti- 
tion, which stabilizes the bubble against Kelvin-Helmhoitz 



and Ray leigh- Taylor instability ( |Dursi|20"07 Dursi &; Pfrom 
~||2008| ). Secondly, if the AGN jet which blew the bubble 



carries large-scale magnetic field, this will stabilize the bub- 
ble. 

The presence of large scale magnetic fields in jets is 
a natural consequence of magnetic launching. It is gener- 
ally accepted that AGN jets are produced (accelerated and 
collimated) by large scale electromagnetic forces originating 
either near the central black hole {e.g. Blandford Sz Znajek| 
1977) or above the accretion disk ( |Blandford &; Payne|1982| 



Lovelace et al.|1987 ). It is expected that after the end of the 



active jet phase, the magnetic field in the cavity will settle 
down to a stable equilibrium. In this paper we discuss the 
properties of such a state. Toroidal fields are dominant in 
jets after expansion takes place, however a purely poloidal 
or toroidal field is unstable; therefore it will relax to a struc- 
ture which contains both components. Thus, the problem of 
ICM bubbles can be split into two parts, the first one is the 
structure of a stable cavity; the second one is the behaviour 
of this cavity as it rises. In this paper we are studying the 
structure of a static cavity after the relaxation has taken 
place. 

As a model problem we are looking for an axisym- 
metric structure of magnetic field within a spherical cav- 
ity. Axisymmetric equilibria are not only simpler to con- 
struct analytically, but there is also good reason to believe 
that they are the simplest and most fundamental configura- 



2 



tions, as the effect of rotation leads to an axially symmet- 
ric sturcture (Lynden-Bell 2006). In the case of fluid stars, 
very general arguments prove that both purely poloidal and 
purely toroidal magnetic fields are unstable (Tayler 



1973 



Flowers & Ruderman 1977). As a result, stable equilibria 



must have considerable toroidal magnetic field component 
(c/. Prendergast 1956). Using numerical methods, the re- 
laxation of an arbitrary magnetic field in a star into an 



axisymmetric equilibrium has been studied (Braithwaite & 
Nordlund||2006| |Braithwaite||2009| see also |Lyutikov||2010| ), 



and recently Braithwaite (2010) has used similar methods to 



study the formation of equilibria in X-ray cavities, finding 
that an arbitrary magnetic field inside a bubble can relax 
into an axisymmetric equilibrium with poloidal and toroidal 
components of roughly equal strength. A rising bubble shall 
experience a stronger stagnation pressure at the top and 
a weaker pressure at the side leading to deformation and 
probably flattening of the bubble. 

Here, in particular, we are looking for solutions with 
vanishing total magnetic field on the surface of the bub- 
ble. This requirement eliminates unphysical surface currents. 
This excludes, in particular, analytical force-free solutions, 
like the spheromak solution. Spheromaks are solutions of 



force-free magnetic fields in axial symmetry ( Chandrasekhar 
&; Kendall 1957]). They have been invoked before in the con- 
text of magnetic bubbles ( jTang||2008| ). However, these so- 
lutions have one particular drawback that excludes them 
as reasonable solution for magnetic field structure in ICM 
cavities: they cannot have both poloidal components equal 
to zero on the surface. Thus, connecting them to unmag- 
netized plasma outside requires surface currents. For exam- 
ple, for the basic spheromak solution the surface current is 
oc sin#. In addition to that, the pressure on the surface of a 
spheromak varies from zero on the axis and goes to a max- 
imum value on the equator, embedding this structure in a 
medium with constant pressure will lead to its deformation. 
It is possible that the system achieves equilibrium after the 
deformation, this problem has been treated numerically by 



Braithwaite (2010) 



Since spheromak solutions correspond to linear force- 
free fields, a large number of high multipole solutions can 
be chosen to construct a system with nearly isotropic sur- 
face pressure. Intracluster cavities are blown by AGN jets 
which may carry large scale magnetic fields, nearly all mod- 
els of AGN jet launching and acceleration require the pres- 
ence of large scale magnetic fields near the central black hole 
( Blandford fe Znajek|1977||Gammie et al.|2003| ); we assume 
that large scale fields survive up to the large distances. The 
jet carries mostly toroidal magnetic field; after the ending 
of an intermittent AGN activity, the magnetic field struc- 
ture within the bubble should reach a stationary state ap- 
proximately on one Alfven crossing time. We assume that 
during this relaxation the magnetic field within the bubble 
preserves a large scale component while small scale fields 
will behave similarly to an additional pressure component. 
In this paper we discuss the structure of magnetic field in 
these relaxed bubbles. 

The magnetic energy in the bubble must be confined by 
the gas, whose pressure must therefore be lower inside the 
bubble than outside. Either the Lorentz force and pressure 
gradient force balance each other throughout the volume of 
the bubble, or the forces act only at the surface of the bubble 



where a singularity in the Lorentz force (a current sheet) is 
balanced by a discontinuity in the gas pressure, in which case 
the interior of the bubble can be force- free. Whilst the latter 
is necessary in the case where the magnetic field dominates 
over the gas pressure in the interior, i.e. as plasma-/^ — »• 0, it 
seems physically unlikely in the high-/3 case - therefore we 
consider here non-force free solutions with no current sheet 
at the boundary. The confinement in this type of solutions 
is not localized on the surface of the system as it has been 
discussed for the case of jets ( | Spruit 1 2 010) but it is extended 
in the whole volume of the system. 

Other options we have attempted, like a purely poloidal 
field with some pressure, in addition to the stability is- 
sues, require surface currents. Our solution has a magnetic 
field with both poloidal and toroidal components and also a 
plasma pressure. We describe the solution and its physical 
meaning in the following section. 



2 MAGNETIC FIELD IN THE CAVITY 

2.1 The Grad-Shafranov equation 

In MHD equilibria the Lorentz force is balanced by pressure 
gradient; we shall assume that gravity is not important as 
we constrain our interest to magnetic cavities and not to 
self- gravitating systems such as magnetic stars, so the equi- 
librium condition is 

JxB = Vp. (1) 
Axially symmetric magnetic fields can be written as 



(2) 

/ Sill V 

where 2ttP is the poloidal magnetic flux and cl is the 
poloidal current enclosed by an axially symmetric loop (e.g. 
Shafranov 1966); we set speed of light to unity. Stationary 
solutions require I — /(P), and the force balance |l]) gives 
the Grad-Shafranov equation 



A*P + F{P)r 2 sin 2 + G{P) = 
where 



A* 



Or 2 



86 



(3) 



(4) 



is the Grad-Shafranov operator, G(P) = All' , I — I(P) and 
F — 47rdp/dP are functions of P. 



2.2 Solution 

We are looking for solutions of (|3| with vanishing magnetic 
field on the surface, as we have assumed the that outside 
plasma is unmagnetized. This requirement leads to an un- 
usual mathematical problem: we seek solutions of an ellip- 
tical equation Q of a function P while matching both its 
value and its derivative on the boundary. 

In other words, on the boundary both the flux function 
P and its normal derivative d r P should be continuous. This 
makes an elliptical equation (J3|, formally, over- constrained, 
as the solution should satisfy both Neumann and Dirichlet 
boundary conditions. 

On the other hand, the Grad-Shafranov equation ([3]) 
has two unknown function F and G, which, as we show 
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below, can be chosen in such a way to satisfy the over- 
constraining boundary condition. This calls for an unusual 
mathematical problem, to find the equation of the type J3|, 
where F(P) and G(P) are functions of the solution only, 
so that the solution itself has to satisfy both Neumann and 
Dirichlet boundary conditions. 

We have found that a linear dependence of the pressure 
and toroidal field function / on P, lead to solutions of the 
desired properties, thus we choose F(P) = Fo and G(P) = 
a 2 P. The angular part of the differential operator of the 
Grad-Shafranov equation admits eigenfunctions which are 
expressed in terms of Legendre polynomials in the following 
way 



Ei sin 



dP;(cos6>) 
d(cos<9) 



(5) 



The fact that the pressure term is multiplied by sin 2 6 con- 
strains our solutions to dipole solutions, any higher multi- 
pole leads to imbalance between the Lorentz and the pres- 
sure forces, for instance see appendix A of Gourgouliat os fc"| 
Vlahakis (2010). Following the method of separation of vari- 



ables for P — sin 2 0f(r) the sin terms drop out and we are 
left with an ordinary differential equation for f(r) 



f" - ^f + F r 2 +a 2 f - 



0. 



(6) 



This equation admits analytical solutions 
sin(ar) 



/ sm ar \ 
j = c\ I a cos(ar) — -J 



> (a 



+C2 a sin(ar) + 



cos(ar) 



r 
F 



(7) 



By demanding that the flux is not infinite at the origin we 
set C2 = 0, and then we observe that the solution we have is 
a spheromak superposed with a uniformly twisted magnetic 
field. 



/ , , sin(ar)\ 
ci I acos(ar) J 



Fq 2 
— r 

a 2 



(8) 



The next step is to enforce the boundary conditions. The 
flux is confined within a sphere of unit radius ro = 1, this 
translates to /(l) = 0. We determine Fq 



Fq = cia 2 (acosa — sin a) 



(9) 



Then we impose the condition that there are no surface cur- 
rents, which is /'(l) — 0, so a is given by 

3a 



tana = 



(10) 



The smallest positive root we find is a — 5.76 and it is 
independent of the choice of normalization, for this root the 
f(r) is zero at r — and r — 1 and it is non-zero in the 
interval between, see Fig. 1; had we chosen a greater root 
for a we would have taken a solution where the flux becomes 
zero within the cavity as well being zero at r — and r — 1. 
For this root of a we find that Fq — 182.67ci. The explicit 
form for the field is the following: 



B(h 



2 cos# . 
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r 

a sin 
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Figure 1. Plot of the solution as it appears in expression |8j 
(solid line), the spheromak component of the field (dotted line) 
and the twisted field (dashed line). 
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Figure 2. Plot of the field components B r , Bq and B^ along r 
for 6 = (the axis), solid line, = 7r/4, dashed line and 6 = tt/2 
(the equator). All three components go to zero at r = 1 and in 
addition B r and B^ reach r = 1 with zero derivatives. 



where the function f(r) and a are determined as we have 
described above, Fig. |2|. 



3 PHYSICAL IMPLICATIONS 

We shall examine the relative contribution of the forces. The 
various terms in the Grad-Shafranov equation ^ are pro- 
portional to the forces acting in the system. The term arising 
from the action of the differential operator A* is due to the 
poloidal field force, the one containing F(P) is due to the 
pressure gradient and finally the one containing G(P) is due 
to the toroidal field force, Fig. (|3j. 



-15 - 

Figure 3. Plot of the relative intensity of the terms that appear 
in the Grad-Shafranov equation |3| in a section of the system 
at = 7r/4. The solid line corresponds to term arising from the 
pressure, the dotted line to the poloidal tern and the dashed line 
to the toroidal term. The three terms add up to zero. 



The solution we have found determines only the pres- 
sure gradient and not its actual value. The pressure profile 
of the system is related to the solution through the relation 
F = 47rdp/dP, given that we know F the actual pressure 
is p = po + FqP/ (47r), po is a constant corresponding to the 
ambient pressure of the medium surrounding the cavity. The 
term FqP is always negative, irrespective of the normaliza- 
tion, thus the pressure always has a dip inside the cavity 
and it requires a positive term po ^ (— FoP/(47r)) m ax, so 
that we never have a negative energy density. The gas pres- 
sure is plotted in Fig. Q. 

We have evaluated the energy of the system as a func- 
tion of a, while keeping the radius fixed. Our aim is to find if 
the case without surface currents has any particular signifi- 
cance in terms of energy. The energy comprises three terms, 
one due to the magnetic field, a second one that is related 
to the term Fo and a third term which is due to the ambient 
pressure po- Although the force equation and the solution 
do not contain the polytropic index (7), this index appears 
in the energy calculation. We are evaluating the quantity W 
which is the energy difference from the energy of a sphere 
of the same radius containing only gas of pressure po 

W = i- J B 2 dV + J (p - Po )dV . (14) 

A case of particular interest is for 7 = 4/3, as the magnetic 
field behaves like a fluid of 7 = 4/3. We have evaluated 
W subject to given helicity, for a range of values of a, but 
keeping the form of Fo so that the system is confined within 
a sphere of given radius, see Fig. 5. For particular values of 
a the surface current vanishes. 

We have found that when 7 = 4/3 the case without 
surface currents is the one that corresponds to minimum 
W which is actually equal to zero, meaning that the energy 
contribution of the magnetic field and of the term containing 
Fo are equal and opposite. This is a consequence of the fact 
that the magnetic field has an effective 7 = 4/3. Had we 
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Figure 4. The gas pressure p as function of r; the solid line is 
for = 7r/2, the dashed line is for = tt/4 and the dotted line is 
the pressure along the axis for 6 = 0. We have chosen po so that 
the minimum of the pressure is zero. 




a 



Figure 5. W as a function of a for a sphere of given radius, flux 
and helicity; and for a fluid of 7 = 4/3. The units are arbitrary. 
The minima occur at a which lead to systems without surface 
currents. The first minimum is at a = 5.76 which is the case we 
study in more detail. 



chosen an other polytropic index for the fluid the partition of 
energy in the fluid and the magnetic field would be different. 

In general, pressure balance requires that any equilib- 
rium satisfies 

^ J B 2 dV + y (P - Po)dV = 0, (15) 

since, averaged over all directions, the magnetic field exerts 
a pressure equal to one third of its energy density - just like 
any other relativist ic fluid. This is confirmed by numerical 
experiments 
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4 SIMULATION 

Numerical simulation is the most straightforward method of 
testing the stability of this magnetic field configuration. In 
this section, such simulations are described. 



The code used is the STAGGER CODE (Nordlund & 
Galsgaard||1995| |Gudiksen fe Nordlund||2005| ), a high-order 
finite- difference Cartesian MHD code which uses a 'hyper- 
diffusion' scheme, a system whereby diffusivities are scaled 
with the length scales present so that badly resolved struc- 
ture near the Nyquist spatial frequency is damped whilst 
preserving well-resolved structure on longer length scales. 
This, and the high-order spatial interpolation and deriva- 
tives (sixth order) and time-stepping (third order) increase 
efficiency by giving a low effective diffusivity at modest reso- 
lution (144 3 here). The value of the diffusivity is set to some- 
what more than the minimum required to avoid numerical 
instability in a variety of situations; since the scheme is de- 
signed primarily to damp behaviour at or near the Nyquist 
frequency, its effect on global structures is minimal - the 
magnetic diffusivity present has a negligible effect on the 
magnetic energy during the simulations. In effect, the dif- 
fusive timescale in this static configuration is at least three 
orders of magnitude greater than the Alfven timescale. The 
code includes Ohmic as well as thermal and kinetic diffusion. 
The use of Cartesian coordinates prevents any problem of 
singularities and simplifies the boundary conditions: peri- 
odic boundaries are used here. 

The numerical model is set up as follows. The compu- 
tational box is a cube of side 3R where R is the radius of the 
bubble. The boundaries are sufficiently far from the bubble 
that they have no significant effect. The bubble is embed- 
ded in an ambient pressure po and can have arbitrary field 
strength up to some maximum where the plasma- /3 is of or- 
der unity; for convenience we choose a field strength such 
that the minimum gas pressure inside the bubble is half 
that outside, i.e. p m i n = po/2. So that all potentially unsta- 
ble modes are excited, a random perturbation of strength 
1% containing a range of all relevant length scales (with a 
flat spectrum) is given to the pressure field. The fluid is then 
evolved forwards in time: stability can be assumed if noth- 
ing happens during a number of Alfven timescales, since any 
instability should grow on this timescale. 

The computational setup and magnetic field are illus- 
trated in Fig. [6] There is no appreciable change in the mag- 
netic field during the simulation. In Fig. [7] are plotted the 
time evolution of the kinetic, thermal and magnetic energies. 
Clearly, the kinetic energy does not grow and none of the 
excited modes is unstable. Since all imaginable modes are 
excited initially, the only logical conclusion is that the con- 
figuration described in the previous section is indeed stable. 

As mentioned above, there is a limit on the strength of 
this field; if the field strength is increased at constant ambi- 
ent pressure po then eventually zero gas pressure is reached 
on the neutral line. This happens when the mean plasma- f3 
in the bubble is approximately unity. Configurations with 
low-/3 must have a different form: since the interior of any 
low-/3 bubble must be roughly force-free, the magnetic field 
is confined by a boundary layer with a strong pressure gradi- 
ent. We can investigate this type of equilibrium by perform- 
ing a simulation similar to that above but with reduced gas 
pressure in the bubble interior. The simplest setup is to use 




Figure 6. Three-dimensional rendering of the equilibrium field 
embedded in the computational box. 
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Figure 7. Kinetic, thermal and magnetic energies as a function of 
time. Time is expressed in units of the Alfven crossing time. The 
initial growth in kinetic energy comes from the perturbation to 
the pressure field; the lack of any further growth implies stability. 



for the initial conditions a magnetic field of the same form 
as before but with increased strength and a low, uniform 
gas pressure inside the bubble, whilst satisfying the pres- 
sure balance condition (15), and to follow the relaxation of 
the configuration into ; 



low-/3 equilibrium. In the following 
simulation, the bubble has a gas pressure of O.lpo so that the 
isotropic magnetic pressure B 2 /24tt — 0.9po- Of course, this 
is a rather crude way of producing a low-beta equilibrium 
and a more realistic method would involve a self-consistent 
cooling mechanism; however the details do not affect the 
qualitative conclusions. 

When the simulation is set running, dynamical motions 
appear and an equilibrium is approached. To aid relaxation 
into the new equilibrium, an artificial friction force of ap- 
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propriate magnitude is added to the momentum equation 
so that oscillations about the new equilibrium are damped. 
In Fig. [8] the structures of the original high-/3 and the new 
low-/3 equilibria are compared. The shape of the bubble has 
deviated from spherica(j and the bulk of the interior is now 
approximately force- free. A boundary layer of high pressure 
gradient confines the magnetic pressure, and since the gas 
pressure must be constant along field lines, there is a narrow 
region of high gas pressure around the axis of symmetry. 
Various quantities are plotted in Fig. [9] on the equatorial 
plane. 

Physically it seems unlikely that a \ow-/3 equilibrium 
could form directly from relaxation of the AGN eject a or 
whatever inflated the bubble. This is because during the re- 
laxation from an arbitrary or 'turbulent' injected magnetic 
field, most of the magnetic energy is transferred to the gas, 
going into some combination of thermal energy and relativis- 
tic particles; in any case, a high-/3 equilibrium is the result. 
A \ow-/3 equilibrium could however form later if the bubble 
gas is able to cool radiatively, in fact not unlike the forma- 
tion in this simulation except that a gradual loss of thermal 
energy would cause quasi-static evolution of the equilibrium 
rather than the sudden change effected here. 



5 OBSERVATIONAL PROPERTIES 

In this section we evaluate some observable properties of the 
system. ICM cavites are observed as depressions in the X-ray 
emission and also are associated with radio lobes, in addi- 
tion other simulated works ( |Dong &; Sto ne 2009; Reynolds 
et al. 1 20051 have similar predictions. For that reason we have 
evaluated the rotation measure of the Faraday rotation the 
cavity will cause to some radiation that lies behind it, the 
synchrotron radiation and its polarization created by the 
cavity and the depression it causes to the X-ray emission. 

There are two ways to test the magnetic field struc- 
ture in the cavity. In the case of a cavity filled with non- 
relativistic material, the structure will cause Faraday rota- 
tion to the radiation that comes from behind the cavity. This 
case is investigated in section |5.1| If the cavity is filled with 
relativistic electrons they will produce synchrotron radiation 
which shall be polarized, in addition relativistic particles do 
not contribute to the rotation measure. Finally we present 
synthetic X-ray images of the system. 

The solution we have found, contains information about 
the magnetic field and consequently about the pressure, 
however in order to evaluate the observable quantities we 
need to know the density and the temperature. It is indeed 
true that the properties of the cavities can be explained by 
assuming a much higher temperature and a low density in- 
side the cavity. However this will allow a great freedom in the 
possible combinations which could lead to similar profiles. 
In our work, we have decided to assume that the density and 
the temperature are associated to the pressure both inside 
and outside the cavity by the same law. Therefore, in the 




Figure 8. Upper panel: cross-section of the high-/3 equilibrium. 
Black lines represent poloidal field; on the right-hand-side the blue 
shading represents the toroidal field (multiplied by the cylindrical 
radius r sin 6) and on the left-hand-side the red shading represents 
gas pressure. The thick yellow line represents the boundary of 
the bubble. That both B^r sin 6 and p — po are functions of the 
flux function and are therefore constant on poloidal field lines is 
clear. Lower panel: the \ow-(3 equilibrium. The pressure is very low 
throughout the bulk of the bubble and there is a steep pressure 
gradient near the boundary which balances the strong Lorentz 
force in that region; note that the field strength now no longer 
goes to zero so gently at the boundary. Note also that during 
relaxation the shape of the bubble has changed somewhat to an 
oblate spheroid, although the volume has remained constant. 



1 There is probably a general tendency for axisymmetric equilib- 
ria to be oblate spheroidal, as this shape allows the poloidal lines 
to be more nearly circular, circular being a lower energy state. 
A proper investigation of non-spherical equilibria is left for the 
future. 
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Figure 9. Upper panel: values of the flux function P, toroidal 
field x cylindrical radius B^r sin and pressure difference p — po 
in the equatorial plane of the high-/3 equilibrium. Middle panel: 
the initial conditions of the low-/3 simulation. Lower panel: the 
new \ow-/3 equilibrium after relaxation. 



case of the rotation measure we shall assume that the non- 
relativistic plasma is described by 7 = 5/3, whereas in the 
case of the synchrotron radiation and the synthetic X-ray 
images, the system is described has 7 = 4/3 and tempera- 
ture depends on the pressure. 

An other issue related to the observability of the quanti- 
ties discussed above, is due to synchrotron emission and ro- 
tation measure from the ICM itself. The intracluster medium 
is denser and cooler compared to the cavity, however accord- 
ing to the magnetic description of cavities we are proposing 
in this paper, the magnetic field inside the cavity shall be 



stronger than the magnetic field of the intracluser medium, 
which in our idealised solution is negligible. In addition if 
we integrate the emissivity along the line of sight the typi- 
cal cavity size shall contribute ~ lOkpc whereas the length 
of the intracluster medium the line of sight crosses, is an 
order of magnitude larger or more. Thus we have compet- 
ing factors, as the stronger magnetic field shall enhance the 
cavity contribution, whereas the lower density and shorter 
length shall depress it. 



5.1 Rotation measure 

We have evaluated the rotation measure corresponding to 
the solutions we have found. The rotation measure is the 
integral of the magnetic field parallel to the line of sight 
multiplied by the electron density along the line of sight. 
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Then we use this information and we draw a map. In this 
problem the magnetic field is described accurately by the so- 
lution. However, in order to estimate the density we need to 
make an assumption for a relation between the pressure and 
the density, in this paradigm we have chosen the density to 
obey an adiabatic relation of 7 = 5/3. We calculate the con- 
tribution to the rotation measure only from the cavity and 
not from the ambient medium. In reality the rotation mea- 
sure may be dominated by the contribution of the ambient 
medium which is more extended and denser. Our prediction 
shall be superimposed to the overall profile of the ICM. 

In addition, as it was described above the pressure is 
determined up to an additive constant, we have chosen this 
additive constant so that the pressure at its minimum inside 
the cavity to be zero, i.e. the f3 « 1 case we have also tested 
cases where the dip in pressure inside the cavity is a fraction 
of the external pressure and we have found that the rotation 
measure is very similar to that of a spheromak configuration 
embedded in a uniform density and pressure environment. 

We have studied three viewing angles with respect to 
the axis of the magnetic field. In the first case we have an 
equatorial view, the axis is perpendicular to the line of sight; 
in the second one the axis is tilted, the angle between the 
line of sight and the cavity axis is 7r/4; and in the third 
case, the polar view, the axis points towards the observer, 



Figs. (10) and (11) 



Coming to the actual value of the rotation measure it 
depends on the density of the electrons, the size of the cav- 
ity and the magnetic field. Setting the maximum value for 
density inside the cavity equal to 10 -4 cra -3 , the magnetic 
field 10/jG and assuming a diameter of lOkpc we find that 
the maximum absolute contribution of the rotation mea- 
sure is about 10 ~ 3 r ad/ cm 2 . Thus if there are any electrons 
available to cause faraday rotation observations of sufficient 
resolution shall see a change in the rotation measure by 



2 x 10 rad/cm . Observations of the Perseus cluster (de 



Bruyn & Brentjens 2005) suggest that these rotation mea- 
sure values are within the limits of observation. 



Figure 10. The rotation measure for three models of cavities 
for two viewing angles, the first row is for the case where the 
line of sight is perpendicular to the axis of the system, whereas 
the second is for systems where the line of sight forms an angle 
of 7r/4 with the axis. The first column corresponds to the case 
where the minimum of pressure in the system is equal to zero. 
The second column is the system corresponds to the case where 
the minimum of the pressure of the system at the minimum is 
half of the ambient pressure. The third column is a spheromak 
filled with gas of constant pressure. 




Figure 11. The rotation measure for the case where the line of 
sight is along the axis of the system. In this case the rotation 
measure is radially symmetric, thus we plot a slice. The squares 
correspond to a solution with ambient pressure such as to go to 
zero at the minimum inside the cavity. The asterisks correspond to 
a solution with an ambient pressure twice as much as the previous 
case. Finally the solid line is the rotation measure for a spheromak 
embedded in a constant pressure environment. 



5.2 Synchrotron emission and polarization 

If the cavity contains relativist ic electrons, the presence of 
the magnetic field will lead to a polarized synchrotron emis- 
sion from the cavity. The intensity of the synchrotron radi- 
ation depends both on the magnetic field and the density of 
the electrons. We assume that the density of the electrons 
is related to the pressure we have evaluated by a relation of 
7 = 4/3, as they are relativistic. The magnetic field is well 
defined by the solution we have found. Again we face the 



Figure 12. First Row: The intensity of the synchrotron emission 
of the system for a system where the line of sight and the axis 
of the cavity form an angle of 7r/2, 7r/4 and from left to right 
respectively. Brighter areas have greater emissivity. Second Row: 
The polarization of the synchrotron emission again for angles of 
7r/2, 7r/4 and 0, the length of the lines is proportional to the 
polarization of the radiation. 



same uncertainties about the minimum pressure in the cav- 
ity, which we now consider to be equal to zero. We find that 
this structure produces synchrotron radiation. The apparent 
profile of the synchrotron radiation depends on the orienta- 



tion of the cavity with respect to the observer, Fig. (12). 

We remark that there are observations which correlate 
low frequency radio emission with x-ray cavities Schwartz 
(2005). The intensity of synchrotron emission (see for in- 
stance Lyutikov et al. (2005}), is given by 



^ + 7/3 / 



K e B 



(n+l)/2 



dl 



(17) 



where n is the spectral index of the population of rela- 
tivistic electrons, in our case we choose it to be n = 2, 
/ is to be integrated along the line of sight, B± is the 
component of the magnetic field which lies perpendicu- 
lar to the line of sight. The function k(v) is defined as 



Kiy) = ^r(22fi)r( 

e is the electron charge 



3n+7 > 
12 ) 



V 27rn 



2 , where 



m e is the rest mass of the electron, 
v is the frequency of the radio emission, and F denotes the 
Euler gamma function. K e can be found from this equation 
n e = K e e~ p dedQ, where n e is the number density of the 
electrons having energy of e to e + de and lie within solid 
angle dQ. Assuming that the maximum magnetic field in 
the cavity is 10/iC, its diameter is lOkpc, its angular size is 
0.5arcmin we find that the maximum intensity of the syn- 
chrotron emission at 100MHz shall be of the order of ~ 5 J. 
This observations shall require resolution of the order of a 
few arcseconds in low frequencies, which is anticipated to be 
within the capabilities of LOFAR jRottge ring 2003). 



5.3 Synthetic X-ray images 

Observations of AGNs have revealed buoyant bubbles as de- 
pressions in the X-ray surface brightness. Following Dong & 
Stone ( 2009 ) we evaluate the X-ray profile of the cavities. We 



assume that the X-ray emissivity of such a system is propor- 
tional to E x r>j p 2 T 1//2 , applying an adiabatic relation with 
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Figure 13. Synthetic X-ray images for a magnetic cavity, the 
line of sight is normal to the axis in the first image, they form an 
angle of 7r/4 in the second one and they are parallel in the third. 
The white areas have lower X-ray brightness. The overall shape 
of the X-ray depression varies as the orientation changes, from a 
system with two lobes to an axially symmetric ring structure. 



an index 7 = 4/3 as we have done in the rotation measure 



we find E x 



(3+ 7 )/(2 7 ) 



Then we integrate E x along the 



line of sight, taking into account the fact that the pressure in 
the external medium is constant and equal to po . Following 
this process we have constructed synthetic X-ray images for 
three different orientations so that the axis and line of sight 



have an angle of 7r/2, tt/4 and 0, Fig. (13). We remark that 



the shapes of the cavities vary from elliptical to spherical 
depending on the orientation, although their boundaries are 
always spherical. 



6 DISCUSSION 

The X-ray cavities appear as areas of lower X-ray emission in 
the intracluster medium, and they originate from AGN jets. 
Purely hydrodynamical simulations in general predict ter- 
minating shocks which are not observed in these systems, in 
addition they are vulnerable to instabilities and demonstrate 
ripples near their edges (|Reynolds et al.|2002). MHD simu- 



lations (Naka mura et al.|2006 ) using the context of magnetic 
towers (Lynden-Bell 2003) suggest that the expansion of a 
magnetic jet in a background density and pressure predict 
the formation of collimated jet and a lobe in the top. |Dong| 
|fe Stone] ( |2009| ) have investigated models of buoyant bubbles 
for a variety of magnetic field configurations and viscosities, 
their simulations suggest that a toroidal field confined in- 
side the cavity is the most promising, however the overall 
structure of the cavity changes drastically as it rises. 

In our study we do not discuss the process of the ini- 
tial formation of the cavity but we focus on the equilibrium 
state after the inflation has taken place. A plausible phys- 
ical description is the following: magnetic jets, after they 
have expanded, are dominated by toroidal field, which is 
now unstable. We suggest that through these instabilities 
the topology of the field changes and the cavities containing 
both toroidal and poloidal field are created. 

Then the cavities rise and expand in the intracluster 
medium, and there is a relation between the radius of the 
cavity and distance from the origin which is discussed in 
detail by Diehl et al. (2008). Our work describes a static 
cavity, a rising cavity will encounter a varying background 
pressure, however this does not pose an issue for our work, 
as it is feasible to balance it with an external environment 
by choosing appropriate values for a and Fq, provided that 
the expansion takes place slowly enough so that the system 
passes through states of equilibrium. However, a possible 



issue could be a pressure environment where the difference in 
background pressure is large enough over lengths of the size 
of the cavity, this may lead to deformation of the structure. 

In our model we have prevented the formation of surface 
currents, this was done to achieve a constant pressure on the 
surface. However if we had allowed the formation of surface 
currents following a force- free model, the cavity would have 
deformed, because the currents vary with the polar angle, 
thus the pressure in the inside of the cavity is not constant 
and cannot balance with a constant external pressure. 

In addition we remark the difficulty of the simulations to 
treat surface currents, as they formally are a discontinuity 
which leads to tedious calculations. Observations are not 
conclusive whether there are surface currents or not, however 
the presence of surface currents has specific disadvantages 
we have discussed previously. We suggest that his solution 
can serve as an initial condition for simulations of rising 
bubbles. 

Coming to the issue of observable properties, they are 
consistent with the structures observed in X rays. However, 
X-ray observations do not give a conclusive answer on the 
presence of magnetic fields. Radio observations of polarized 
radio emission or rotation measure of the system can give a 
more clear answer the question of the magnetic field. How- 
ever, we are not aware of such measurement of sufficient res- 
olution to determine the radio behaviour of the bubble. We 
remark however that magnetic fields are inferred in large- 
scale structures coinciding with the Hercules and Perseus- 



Pisces superclusters (Xu et al. 2006) and also a study of 



the polarization of Perseus cluster by 



(2005) 



de Bruyn Brent jens| 



We remark that similar structures occur as eject a in 
magnetar giant flares (Gae nsler et al.|[2005 ). They are cre- 
ated by strong magnetic field and plasma which is driven by 
the magnetic action on the magnetar, but they are scaled 
to smaller sizes. As our solution does not have an intrinsic 
scale can be applied to those configurations. 



7 CONCLUSIONS 

In this paper we have found analytical solutions for mag- 
netic cavities without surface currents. The cavities contain 
a magnetic field with poloidal and toroidal components and 
a hot plasma. The structure of the fields is such so that 
they drop gradually to zero at the end of the cavity. Because 
of that, there are no surface currents, unlike the force- free 
fields which require surface currents. There is a strong qual- 
itative similarity between this field configuration and that 



found with numerical methods by Braithwaite (2010). The 



non- force- free equilibrium found here must of course have 
a plasma- /3 higher than about unity; a \ow-/3 equilibrium 
must be force-free in the interior with a current sheet at the 
boundary. 

A numerical test has confirmed the stability of the sim- 
plest (fundamental radial mode) equilibrium. Furthermore, 
simulations have been used to examine the behavior of the 
equilibrium when the gas pressure drops by a large factor 
(via cooling) and takes the bubble from the f3 y 1 regime 
into the f$ ■< 1 regime. It is found that, as expected, the 
bubble becomes approximately force-free in the bulk with a 
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thin layer of high current density at the boundary and along 
the axis of symmetry. 

These structures correspond to rotation measures which 
are distinguishable from other structures and especially 
spheromaks. These are clear when the line of sight is nor- 
mal to their axis, however they get saturated if the if the 
magnetic field is too weak. When the angle between the 
line of sight and the axis of the system is small they resem- 
ble spheromak structures more closely, but there are distin- 
guishable properties visible with a sufficient resolution. In 
addition we have constructed synthetic X-ray images that 
show remarkable resemblance to x-ray cavities observed. 

We remark that although a bubble has a spherical 
boundary it may appear as an ellipsoid in the X-ray ob- 
servations, depending on the orientation. This accounts for 
the variety of the observed shapes. In any case the cavities 
may have intrinsically various shapes. In addition we have 
evaluated the rotation measure that will affect a background 
source if they contain non-relativistic plasma and the syn- 
chrotron radiation and its polarization if they contain rela- 
tivists particles. 
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